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We realize tunable coupling between two superconducting transmission line resonators. The cou¬ 
pling is mediated by a non-hysteretic rf SQUID acting as a flux-tunable mutual inductance between 
the resonators. From the mode distance observed in spectroscopy experiments, we derive a cou¬ 
pling strength g/27T ranging between —320 MHz and 37 MHz. In the case of ^~0 the microwave 
power cross transmission between the two resonators can be reduced by almost four orders of mag¬ 
nitude compared to the case where the coupling is switched on. In addition, we observe parametric 
amplification by applying a suitable additional drive tone. 


I. INTRODUCTION 


In circuit quantum electrodynamics, the controllable in¬ 
teraction of circuit elements is a highly desirable resource 
for quantum computation and quantum simulation ex¬ 
periments. The most common method is a static capaci¬ 
tive or inductive coupling between cavities and/or qubits. 
In such a system, photon exchange can be controlled by 
either tuning the circuit elements in and out of resonance 
or using sideband transitions.While this approach has 
proven to be useful for few coupled circuit elements, it 
seems impracticable for larger systems, where it is hard to 
provide sufficient detunings between all circuit elements.^ 
Therefore, one may alternatively use tunable coupling 
elements such as qubits^“^ or SQUIDs.^^“^^ One partic¬ 
ular example for an interesting application of such ac¬ 
tively coupled circuit elements are quantum simulations 
of bosonic many-body Hamiltonians.In such a sce¬ 
nario, the bosonic degrees of freedom can be represented 
by networks of (possibly nonlinear) superconducting res¬ 
onators. For this quantum simulator, a tunable coupler 
would constitute an important control knob. A more 
general scope of this device is the controllable routing of 
photonic states on a chip, which is interesting for quan¬ 
tum information as well as quantum simulation experi¬ 
ments. 

In this work, we experimentally investigate the case of 
two nearly frequency-degenerate superconducting trans¬ 
mission line resonators coupled by an rf SQUID in the 
spirit of Refs. [20 and 21]. The rf SQUID, a supercon¬ 
ducting loop intersected by a single Josephson junction, 
acts as a flux-tunable mutual inductance, mediating a 
flux tunable coupling between the two resonators. Tun¬ 
ing the rf SQUID by an external flux, we can vary the 


coupling strength between —320 MHz and 37 MHz and 
observe the corresponding coupled modes in our system. 
Furthermore, we find that the coupling can be suppressed 
by approximately four orders of magnitude and, there¬ 
fore, efficiently be switched off. Finally, the flux-tunable 
inductance provided by the rf SQUID gives access to 
parametric amplification with gains of up to 20 dB. Our 
results are particularly relevant for the realization of scal¬ 
able arrays of superconducting resonators for quantum 
information processing and quantum simulation applica¬ 
tions. In this context, we emphasize that the rf SQUID 
coupler requires only a single Josephson junction and a 
single control field per coupler. 


II. SYSTEM HAMILTONIAN 

An optical micrograph of the sample is shown in Fig. 1(a). 
Our system is comprised of an rf SQUID galvanically 
coupled to the center conductor of two coplanar stripline 
resonators. These resonators, A and B, can be described 
as quantum harmonic oscillators using the Hamiltonian 

iJres = huJAd^Ci + . (1) 

Here, uja and cjb are the resonance frequencies and a^, , 

a, and b are the bosonic creation and annihilation opera¬ 
tors. The effect of the rf SQUID on the system properties 
can be modeled in terms of an effective inductance. 

The fluxes and <Fb generated by the resonators in 
the rf SQUID give rise to an inductive interaction en¬ 
ergy. The rf SQUID consists of a superconducting loop 
with inductance Lg, which is interrupted by a Josephson 
junction with critical current Jc- The flux ^ threading 
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FIG. 1. Sample and measurement setup, (a) Sketch of the 
sample chip. Dark blue: Resonator groundplanes. Green: 
Resonator center conductors; light blue: Feed line center con¬ 
ductors. The insets are (false color) optical micrographs of 
the coupling capacitors (orange), the rf SQUID (red), and 
the rf SQUID junction (purple), (b) Sketch of the device mea¬ 
surement setup indicating the attenuated input lines (wiggly 
arrows), the output lines including cryogenic and room tem¬ 
perature microwave amplihers (triangular symbols), and pos¬ 
sible measurement paths (blue dashed arrows). 


the SQUID loop gives rise to a circulating current 


Calculating the injected fluxes in Eq. (5) as described in 
Ref. 25, two essential properties of the system become 
obvious. First, the g-terms in Eq. (5) lead to dressed 
resonator frequencies 


<^A,B = ^A,B 


~ ^A,B 


Q -2 


1 - 




^A,B^rf(^) 

Ll 


^A,B^rf(^) 


(6) 


where I/a,b is the inductance of the resonators and Lc 
the inductance of the segment shared between resonator 
and rf SQUID. The second effect of Eq. (5), caused by 
term oc^a^b, is a flux dependent coupling 


9*b(4>) = 

between the resonators. Due to their vicinity on the chip, 
there is a also flux independent direct inductive coupling 
component go between the resonators. Thus the total 
coupling reads 


^(^) — ^ab(^) + ^0- (8) 

Equation (4) shows that ^ab(^) can be positive or neg¬ 
ative depending on the applied flux. By applying a suit¬ 
able flux, the rf SQUID mediated coupling compensates 
the direct inductive coupling. In this way one can turn 
on and off the net coupling between the resonators. Af¬ 
ter a rotating wave approximation the full Hamiltonian 
reads 


/s($) = -/cSin(27r$/$o), (2) 

where 4>o is the flux quantum. Here, 4> is the sum of the 
externally applied flux 4>ext and the flux generated by A, 

i> = i>ext+is/s(^). (3) 

Since, in the experiment, the screening parameter 
(3 = 27rI/s/c/4>o < 1, the dependence of the total flux on 
the external flux is single-valued. Erom Eq. (2) and 
Eq. (3), an expression for the effective SQUID inductance 
is obtained, 

1 dh 1 ^cos(27r^) 

Arf(4>) 94>ext Eg 1 + ^cos (27r^) 

With this result, we obtain the inductive contribution to 
the system Hamiltonian, 

^ _ ($A - ^b)' _ ^a + - 2^a^b 

2Lrf(^) 2Lrf(^) • 



The eigenvalues of Eq. (9) correspond to the new eigen- 
frequencies 

0., = ± . (10) 

III. SAMPLE AND MEASUREMENT SETUP 

Eigure 1(a) shows the layout of the sample chip. In 
the resonator design, we omit the second groundplane 
to reduce the direct geometric coupling between the two 
resonators. The rf SQUID is galvanically connected to 
both center strips of the resonators over a length of 
200 pm. The sample is fabricated as follows. Eirst, a 
100 nm thick niobium layer is sputter deposited onto a 
250 pm thick, thermally oxidized silicon wafer. The res¬ 
onators and the SQUID loop are patterned using opti¬ 
cal lithography and reactive ion etching. The Josephson 
junction of the SQUID is fabricated in a Nb/AlOx/Nb 
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FIG. 2. Through transmission magnitude (color coded) as a function of probe frequency and externally applied flux for (a) 
resonator A and (b) resonator B. The red arrows mark the flux values, for which transmission vs. frequency cuts are shown 
in Fig. 4. (c) Fit (red line) of Eq. (10) to the extracted center frequencies (crosses). The modes are take from the trough 
measurements, where they are more pronounced. This is especially necessary when the coupling is much smaller then the 
detuning of the resonators. 


trilayer process. The resonators have a characteris¬ 
tic impedance of Zq = 64 and the resonance frequen¬ 
cies u;A/27r = 6.461 GHz and ci;B/27r = 6.482 GHz. The 
SQUID loop has dimensions of 200 pm x 100 pm and a 
screening parameter [3 = 0.934 to maximize the coupling 
according to Eq. (4) while keeping the SQUID mono¬ 
stable. The sample is mounted inside a gold plated cop¬ 
per box, which is attached to the base temperature stage 
of a dilution refrigerator operating at 26 mK. A super¬ 
conducting solenoid attached to the top of the sample 
box is used to generate the external flux applied to the 
rf SQUID. As depicted in Fig. 1(b), one port of each res¬ 
onator is connected to an attenuated input microwave 
line, whereas the remaining ports are connected to out¬ 
put lines containing microwave amplifiers. 


IV. RESONATOR SPECTROSCOPY 

We first extract the properties of the rf SQUID coupler 
from transmission measurements through the individual 
resonators. As indicated in Fig. 1, we call this type of 
measurement a “through-measurement”. In contrast, in 
the “cross-measurements” we inject a signal into one of 
the resonators and probe the output of the other one. In 
these measurements, the applied microwave power corre¬ 
sponds to an average photon number of about one in the 
resonators. In Fig. 2(a) and Fig. 2(b), the through mea¬ 
surements of resonator A and B are shown depending on 
the applied flux ^ext- According to Eq. (4) and Eq. (5), 
the modulation of the resonator modes due to the pres¬ 
ence of the rf SQUID is $o-periodic and symmetric with 
respect to <l>ext = ^o/2. The two modes of Eq. (10) man¬ 
ifest themselves as two resonances in the spectroscopy 
data. As expected, we observe a flux dependent mode 
distance, caused by the flux tunable mutual inductance of 
the rf SQUID. For most flux values, one observes two res¬ 
onance peaks independent of the chosen input and output 
port. However near ^ext = 0.468 and ^ext = 0.532 


only a single peak is present in the through measurements 
and cross transmission is strongly suppressed (see Fig. 3). 
These are the points where the SQUID-mediated cou¬ 
pling compensates the direct inductive coupling result¬ 
ing in a vanishing total coupling and, hence, completely 
decoupled resonators. Note that for g{^) = 0, the Hamil¬ 
tonian of Eq. (9) becomes diagonal and each of the two 
modes reflects the excitation of one of the resonators. In 
Fig. 2(c), the center frequencies of the normal modes 
derived from the data in Fig. 2(a) and Fig. 2(b) are plot¬ 
ted along with a fit using Eq. (10). Upon closer inspec- 
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EIG. 3. Gross-transmission from resonator A to resonator B 
as a function of the applied flux. Near 4>ext = 0.468 4>o and 
^ext = 0.5324>o (white arrows), the resonators decouple and 
signal transmission is blocked. The red arrows mark the flux 
values, for which transmission vs. frequency cuts are shown 
in Eig. 4. 

tion, we find a minimum distance on the order of 20 MHz 
between these modes. This finite gap is caused by a small 
detuning A = cja — = 27r x 21.3 MHz of the resonators. 

We also observe different decay rates of the resonators, 
which we extract from the through measurements of 
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FIG. 4. Through- (red) and cross-measurements (blue), obtained for three different flux values, (a) 4>ext = 4>o/2,|^| ^ | A|: 
strongly coupled regime, only the mode is shown, (b) 4>ext = 0.4684>o/2, |p|~0: decoupling point, (c) 4>ext = 0.4394>o/2, 


both resonators at a decoupling point. Lorentzian 
fits lead to 7A/27r = 3.6 MHz and 7B/27r = 6.1 MHz. 
The fact that resonator A has smaller width at half 
maximum and a and a slightly higher eigenfrequency 
than resonator B, indicates a smaller effective coupling 
capacitance,^^ which we attribute to fabrication or sam¬ 
ple contacting imperfections. Therefore, we assume 
La/Lb ~ 1- Furthermore, we define the fitting parameter 
go = yu;A<^B^c/(V^A^BAs). In this way, the rf SQUID 
coupling reads ^ab = flo/3 cos (27r^)/(1+^ cos (27r^)). 
Fitting the data, we obtain ^o/27r = 29.0 MHz and 
5fo/27r = 20.4 MHz. From the mode distance we find 
5f(4>)/27r ranging between 37MHz and —320 MHz. Be¬ 
yond this value, the lower mode becomes too broad due 
to its steep flux dependence. 

Next, we analyze the properties of our device in the 
coupled and decoupled state in more detail. Because 
of the small detuning of the resonators, the coupled 
modes are not necessarily symmetric and antisymmet¬ 
ric superpositions of the uncoupled modes. This is 
also seen in the spectroscopy of the single resonators 
(see Fig. 2), where the modes have different intensities. 
The mode mixing can be estimated from the eigenvec¬ 
tors of the Hamiltonian in Eq. (9). For 5f/27r = 37MHz 
and g/27T = —320 MHz we obtain the mixing ratios 63:37 
and 52:48, respectively. Hence, in the latter case, our 
sample satisfies the condition | 5 f|::^|A|, where the de¬ 
tuning becomes insignificant. In the decoupled case near 
4>ext = 0.468 4>o and 4>ext = 0.532 4>o, the off-diagonal ele¬ 
ments in the Hamiltonian of Eq. (9) vanish and the modes 
are pure excitations of resonator A or B. 

In this situation, it is particularly instructive to exam¬ 
ine the cross-transmission spectra such as the one shown 
in Eig. 3, where resonator A is driven and resonator B 
is probed. Here, we clearly see that in two narrow re¬ 
gions around 4>ext = 0.4684>o and 4>ext = 0.5324>o, where 
the net coupling g{^) approaches zero, the microwave 
transmission between the resonators is blocked. We gain 
further insight by comparing the through- and cross¬ 
transmission spectra in Eig. 4(a) and Eig. 4(b). Eor 


|^(4>)| ^ |A| [see Eig. 4(a)], both measurements ex¬ 
hibit similar peak heights. Since both measurements 
use the same output line, we relate the small difference 
of approximately 1.5 dB mainly to the slightly different 
losses in the input lines. Eor ^(4>)«0, however, the 
cross-transmission is suppressed by 40 dB on resonance 
as shown in Eig. 4(b), corresponding to a relative trans¬ 
mission change of 38.5 dB. This result confirms that we 
can sufficiently compensate the direct inductive coupling 
with the tunable SQUID-mediated coupling. Einally, in 
Eig. 4(c) we show the transmission for a flux value, where 
^(4>) and A are comparable. In the through measure¬ 
ment, the detuning manifests itself in the form of un¬ 
equal peak heights and an anti-resonance dip, which is 
not centered between the resonance peaks. 


V. PARAMETRIC AMPLIFICATION 

So far, we have treated the rf SQUID as a flux-tunable 
coupler between the two resonators. However, we can 
also interpret it as a flux-tunable inductance, parametri¬ 
cally modifying the eigenfrequency of the modes. Then, a 
harmonic flux drive applied in addition to the static flux 
allows for the study of parametric effects. Our experi¬ 
ments in this direction are inspired by analog experiments 
on Josephson parametric amplifierand converter. 

In order to operate our device as a flux-driven ampli¬ 
fier, we apply an additional drive tone at cjd = 2Ui to 
the input line. The drive tone periodically modulates 
the flux threading the SQUID loop and therefore the 
mode frequency, leading to parametric amplification. To 
characterize the performance of the resulting paramet¬ 
ric amplifier, we calculate the power gain (G) as well as 
the bandwidth characterized by the full width half max¬ 
imum AU of the amplified signals. Eigure 5 shows the 
transmission in the vicinity of Ui/27r = 6.472 GHz and 
4>ext = 0.450 4>o for different values of the drive power. 
While the gain is increasing for higher drive strength, 
the bandwidth decreases as expected. Eor a nonde¬ 
generate (phase-insensitive) gain of G = 20dB, the gain- 
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FIG. 5. Mode for 4>ext = 0.450 4>o with an additional drive 
tone Ct;D/27r= 12.944 GHz for different drive powers Pd- The 
drive power refers to the power at sample input. 


bandwidth product is GAO/27r = 21.5 MHz. Checking 
the theoretical relation^^ AO(G) = const., we find a 
maximum deviation of 2.5 similar to other experiments.^^ 


VI. CONCLUSIONS 

In conclusion, we present a flux-tunable coupling between 
two superconducting resonators based on an rf SQUID. 
Spectroscopically, we measure negative and positive cou¬ 
plings ranging from —320 MHz to 37 MHz. Furthermore, 
the observed suppression of the cross-transmission of up 
to 38.5 dB proves the ability to effectively turn off the 
coupling and is a significant improvement compared to 
previous work.^’^ With this performance, our coupler is 
considered a useful tool both for quantum computation, 
where a controlled nearest neighbor interaction is re¬ 
quired and our coupler could be used to route information 
in a controllable way on a chip and for quantum simula¬ 
tion experiments,where it could be useful not 
only to change the amplitude of the coupling constant but 
also its sign.^^ Finally, we investigate parametric ampli¬ 
fication of the device and observe a nondegenerate gain 
as large as 20 dB. 
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